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ABSTRACT 

The theory of a Gerdien Condenser operating in a collision- 
controlled medium is reviewed. Design and electronics of a Gerdien 
condenser probe suitable for flying on the Areas rocket is presented. 
Aerodynamics properties of the instrument in continuous flow are 
discussed. The method of data reduction and experimental results of 
one successful flight at White Sand Missile Range (WSMR) , New Mexico 
on 11 January 1974 is reported. The result of this investigation 
shows: positive ions in two relatively distinct mobility groups between 
47 and 65 km and a more continuous distribution of mobilities between 


38 and 47 km. 



CHAPTER I 


INTRODUCTION 

1.1 Importance of the Problem 

The ionospheric D-region is located at about 60 to 90 kilometers 
(km) altitude and is largely responsible for the absorption of radio 
waves. Historically, the most important ionospheric measurements in 
this region have utilized remote probing with such waves (Deeks, 1966; 
Coyne and Belrose, 1972; Lee and Ferraro, 1969), However, in terms 
of the physical parameters of this region, these methods provide 
information only about free electrons, specifically density and 
collision frequency . The various remote sounding techniques are 
ineffective below about 65 km; and some of these methods have limita- 
tions in that they can only be used during disturbed conditions, when 
the electron concentration is high enough to give a reliable measurement 
(Thrane et al., 1968). 

Direct measurements in this region may be made only by sounding 
rocket at present, since it is above the region where balloons can 
ascend and below the normal operational altitude of satellites. 

However, the shock waves generally characterized at sounding rockets 
may seriously affect direct measurements in this region (Hale et al., 
1972) . Some parachute borne measurements which avoid this problem have 
been performed, including the blunt probe. 

The basic quantities measured by the blunt probe are the polar 
charged particle conductivities, a+ and C-, and in order to obtain ion 
density (N) the mobility (k) must be known (since a - Nek). Heretofore, 
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positivs ion mobility was Essumsd to bo that spproprifltfi for small 
molecular ions at a temperature and density determined as a function 
of altitude from an appropriate atmospheric model (CIRA) . Recently, 
however; a large amount of evidence has accumulated that these values 
of mobility may not in fact be correct, and that a substantial number 
of ions may not be "small" (AMU <100), but quite large and of very 
much lower mobility; a summary of this evidence has been given by 
Chesworth (1974) . A direct measurement utilizing a subsonic Gerdien 
condenser "aspiration probe" was made by Rose and Widdel (1972) which 
indeed showed positive ions to have such small mobilities, as well as 
mobilities for the "small" ions larger than expected. However, all of 
this data was obtained from a single launch of the payload. 

It can be concluded, therefore; that there is a need for more 
information about ionospheric mobility. The Gerdien condenser is one 
device that can be used to measure charged particle mobility. The 
angle of attack problem associated with earlier instruments (Pederson, 
1964) can be largely eliminated by the use of the more aerodynamically 
stable decelerating devices which are available today. 

1. 2 Previous Related Research 

The Gerdien condenser has been named after Gerdien (1905), who 
first demonstrated use of two concentric cylinders for ion collection 
at ground level. Some investigators have used this instrument to 
measure atmospheric ion conductivities and mobilities for many years. 

A review of their work has been given by Knoll et al. (1964). Further 
development of the Gerdien condenser as a mobility spectrometer have 
been reported (Israel and Schulz, 1933; Nolan and Kenny, 1952). Israel 
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and Schulz were able to differentiate four groups of ions and their 
seasonal and geographical variations by using these devices at sea level. 

’The Gerdien condenser has been used in numerous experiments to 
measure low altitude (up to about 30 km) atmospheric ion density and 
conductivity on balloon flights (Stergis et al., 1955; Woessner et alo, 
1958; Kroening, 1960). The results of these experiments may be subject 
to some serious errors (Paltridge, 1965). Some measurements have also 
been done using the instrument on aircraft (Sagalyn, 1958; Kraakevik, 
1958). Due to aircraft noise there may be some uncertainties in this 
type of measurements. 

Bourdeau, Whipple and Clark (1959) have described a rocket experi- 
ment using two Gerdien condensers mounted on the noset ip of a Viking 
rocket* Their measurements show an ion conductivity less than they 
expected. The authors have pointed out the uncertainty in their 
measurement, since the shock waves may disturb the air flow through 
condenser as the rocket is moving with a supersonic velocity (M = 4) . 

The shock waves may also cause thermal ionization and alter the com- 
position of the medium. Whipple (1964) has described further develop- 
ment of Gerdien condenser systems mounted at the tip of a rocket nose- 
cone, and some of measurement results from Aerobee experiment at Fort 
Churchill, Canada, has been reported by Whipple (1964). 

Pederson (1964) has employed a parachute-borne Gerdien condenser 
technique to measure ionospheric ion densities. In addition to the 
angle of attack problem, the wake of the detached rocket nosecone in 
front of the probe complicates the air flow through the condenser in 
Pederson 1 s experiment. 
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Perhaps thfi most rslxsbl^ lotiosphsxxc chsrs^ pflrtxclft msssur6Qi6iit 
using the parachute-borne Gerdien condenser to date was by Rose and 
Widdel (1972) . In addition to their analyzer system which has a second 
aspiration system behind it to remove all charged particles from the 
air stream (keeping the charge neutrality at the end of the sonde), a 
flowmeter similar (in size and geometry) to the condenser was used to 
determine the velocity of the air stream through the aspiration 
condenser. They have also pointed out that nickel coating of the system 
eliminates a surface effect present with aluminum which caused a zero 
shift in the current-voltage characteristic, and by using care in 
aerodynamical design their system is relatively unaffected to an angle 
of attack up to 20° to 25°. Their experimental results on May 15, 

1968 from Perdasdef ogu, Sardina shows: two groups of positive ions in 
the height region 58 to 72 km, one group of positive ion below about 
52 km, and no data in the altitude range of 52 to 58 km. Their data 
shows only one group of negative ion which have a mobility about the 
same as the mobility of lighter group of positive ions. 

More recently Conley (1974), after a review of Gerdien condenser 
theory, has presented the experimental results of D-region positive ion 
density and mobility measurements during a PCA event. In four experi- 
ments which took place at Churchill Research Range (CRR) , the Gerdien 
condenser was mounted on a black Brant rocket having a supersonic 
speed (M = 3 to 4) . Conley's data shows one group of positive ions 
with mobility slightly less than the lighter ions (between about 58 and 
72 km) and mixed ion (below 52 km), measured by Rose and Widdel (1972). 
However, the data indicate no low mobility group of ions, which may be 
due to the supersonic nature of the experiment; the shock wave may 
prevent the sampling of these ions. 
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CHAPTER II 
THEORY 


2.1 Ion Collection 

The Gerdien condenser consists basically of two concentric 
cylinders. Figure l.a, positioned so that air flows through the 
instrument at sufficiently low potential difference between the inner 
and outer cylinders; ions are collected at a rate proportional to the 
flow rate and this potential difference. 

Before saturation, the general equation for positive ion collec- 
tion can be written as: 

J = -eDVN + + kN + eVV + N + eS (2.1) 

where J is current density; D is coefficient of diffusion of positive 
ion; k is mobility of positive ion, N is number density of positive 
ion; V is potential; e is electrostatic charge; and S is the convection 
velocity. 

The first term in equation (2.1) is due to diffusion. The second 
term represents the contribution of conduction, and the third term is 
due to convection. In the classical operation of the condenser, it is 
usually assumed that air flow is parallel to the axis of the instrument* 
Under this condition the convection term can be neglected and equation 
(2.1) reduces to: 



io di/dvCCNK, 
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Under the assumption that VN + has a value at most of the same 
order of VV, the first term in equation (2.2) can be neglected compared 
to the second term. This can be shown using Einstein’s relation: 


D _ KT 
k e 


(2.3) 


where K is the Boltzmann constant and T is the temperature in degrees 
Kelvin. 

KT 

In upper atmosphere the value of — < -03V, therefore, for a 

e — - 

probe potential much greater than this equation (2.2) becomes simply: 

J = kN + e?V (2.4) 


It has been shown (Hoult, 1965; Sonin, 1967) that the equation of 
the form (2.4) is applicable to a wide class of probes operating in a 
continuous medium, 

2,2 The Current-Voltage Characteristic of the Gerdien Condenser 

Consider the Gerdien condenser shown in Figure l.a consisting of 
two concentric cylinders of radii R q and R „ An electric field is set 
up between the two cylinders in such a manner that the ions to be 
measured drift toward the central electrode which is the collector in 
this arrangement. In order to obtain the current-voltage characteristic 
of the Gerdien condenser, the trajectory of a positive ion, entering 
the condenser at a radius R , would be considered (Pederson, 1964) - 

X 

The radial drift velocity Ur is: 

U r = k + E(r) (2.5) 

where E (r) is the radial electric field and is given by: 
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E (r) = - (2 * 6) 

ln \ 

where r is any radial distance from condenser axis. The volume rate of 

2 2 

air flow through the cross-section tt(R - r ) is: 

X 

R 

R f X 

tj> x = j U x (r)2TTrdr (2.7) 

r 4 

where U (r) is air stream velocity and is a function of x and r 
(Paltridge, 1965). In a time increment dt, the following relation 
holds: 


, dr dx 

dt * r ■ u (O 

r x 


U X (r) 

dx = dr 


substituting the value of from equation (2,5), dx integrates to: 


R 


x 55 


U (r) 
x 

K k + E (r) 


-In 


L 

r 


dr = 


R 


k + V 


— J u x (r)rdr 


( 2 . 8 ) 


R 


x 


x 


combining equations (2.7) and (2.8) 


R 

In — R 

i R n x 

X 21 V r 


(2.9) 


If a positive ion enters the condenser at a radius R xl> a voltage 
is necessary to obtain x = L, at the end of the central electrode, 
when r = R^. All ions entering the condenser at a radius of less than 
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l xl are collected by the central electrode resulting In a current of: 


h - n+ ^r " 2 - N+e f k \ 


( 2 . 10 ) 


OtTC-T, 

where £ is the permittivity of the medium and C = - — is the capacity 

R 

1 -J! 

of condenser* n R 

i 

It can be seen from equation (2*10) that before saturation, the 

31 i 

current is proportion to a voltage V , 3 ^““ proportion to the 
+ + + ^ 

conductivity O = eN k * The saturation current is given by: 


.+ 

1 

o 


1 R 
, T + A O 

= N e C)) R 

i 


N + e £ k + V 
e o 


( 2 . 11 ) 


If instead of U (r) an average velocity U be assumed, then using 
equations ( 2 * 8 ) and ( 2 . 11 ) the saturation voltage can be written as: 

R 




V = 
0 




Ck 


( 2 . 12 ) 


The equations (2.11) and (2.12) are unchanged if the electric 
field is reversed and the outer cylinder becomes the collector. They 
are also valid for negative ions provided that N and k are replaced 
by N and k , respectively. The current voltage characteristic of the 
Gerdien condenser is shown in Figure l.b. One type of positive ion 
and one type of negative ion and electron have been assumed to be 
present in the air aspirated through the condenser. The negative 
saturation current i Q is proportional to U (N~ + N) where N is the 
electron density. Electrons having a high mobility are collected for 
very small voltages as shown in the figure. Stray fields at two ends 
of the condenser have been neglected. 
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The current-voltage characteristics of the Gerdien condenser when 
qjcq different types of xons of different densities and mobilities 
can he obtained by superposition* Under this condition equations 
( 2 . 10 ) and ( 2 . 11 ) become: 

i = e — V 2 k N (2.13) 

z n . n 


i 

o 



n 


UN 
n n 


(2*14) 


Figure l*c shows the ideal current-voltage characteristic of a Gerdien 

condenser 5 if two types of positive ions are present* The saturation 

current i is proportional to (N* + N*)U. The derivative of the current- 

f - j- 

voltage characteristic is proportional to (N^ + N 2 K 2> £or V < V ol 

and proportional to for V ol < V < V o 2 - 

From the foregoing discussion, it is evident that the ion conductiv- 
ity, N- e K— , can be deduced from the slope of .the current-voltage, charac- 
teristic, and U can be derived from the saturation current. In 
order to obtain information about K— either N— or U must be known. The 
Gerdien condenser can be operated as a mobility spectrometer as demon- 
strated in Figure l.c. It is, however, necessary to have an accurate 
current-voltage characteristics if ions of different mobilities are to 
be distinguished. 


2.3 Space Charge Effect 

The relation for electric field in the last section (equation 2.6) 
is a solution of Laplace's equation which assumes the existence of 
space charge neutrality in the Gerdien condenser. If ions of one 
polarity are collected more quickly than the other polarity, a space 
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charge will be present between two electrodes of the Gerdien condenser. 
This space charge may distort the electric field in the condenser. An 
estimate of the space charge can be made by solving the Poisson's 
equation assuming one type of particle remains constant through the 
condenser. Poisson's equation for a Gerdien condenser can be written 
as : 


d 2 V 1 dV _ N + e 

. 2 + r dr e 

dr 


(2.15) 


whose general solution for constant positive charge density is: 

dV 


E (r) = 


dr 


N + er B 
4e r 


(2.16) 


and 


AT + 2 
N er 

V - - 4* B In r + C 

4£ 


(2*17) 


where B and C are the constants of integration and are evaluated by 

applying the boundary conditions at r = R^, V = 0; and at r = R qJ 

V = V : 
o 


B *= 


In 


R 

o 

R , 
x 


2 2 

M + R - R7 
N e o x 

4s R 


In 


R. 

x 


C = 


+ r2 

- B In E. + N e r 
i 4e 


(2.18) 


(2.19) 


Substituting equations (2.18) and (2.19) for B and C, equations (2,16) 
and (2.17) become: 


. V 


E (r) = 


R 


N + e 

2e 


r - R - 


R 2 

i 


R 


( 2 . 20 ) 


r In 


2 r In —■ 
i_ 
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-V 0 1 R. M + , 

V cr 2 

R 4e o 


R 2 , 


In 


R. 

x 


In 


R. 

x 


/ 2 _ 2 n 
(r - R.) 

l 


(2 • 21 ) 


In 


R. 


The distribution of electric field for positive charge densities 
x ^ “3 A -3 "t 

of N = 10 cm , N =10 cm and N =0 with applied potential of 

V ~ 5v and V = lOv is shown in Figure 2, It is shown that for charge 

4 -3 

densities of the order up to 10 cm , the space charge neutrality 

assumption causes no significant error. However, above 80 km, it is 

believed that electron and positive ions are the dominant charge 

particles. Electrons having a high mobility, are picked up by the 

collector at the opening of Gerdien condenser. Positive ions follow 

the air stream through the condenser, where they represent a space 

5 -3 

charge (density 10 cm or more) in the condenser. Therefore, when 
Gerdien condenser is used above 80 km the effect of space charge must 
be taken into the account. 


2.4 Ion Mobility 

When a mixture of ions and gas are under the influence of an 
electric field they gain energy from the field. If the energy acquired 
in one mean free path is small compared to thermal energy: 

(- + eEX « KT (2.22) 

m H 

where M is the mass of the gas molecule 
m is the mass of the ions 
A is the mean free path 

and after substitution inequality (2.22) becomes: 



L3 - 



Figure 


f the Gerdien 
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(— + — ) e E « Pa (2.23) 

V m' 

where P is the gas pressure and O the collision cross section. If 
there is not too much difference between the masses of ions and 
molecules, inequality (2.23) reduces to. 


E 

P 


« 



2 


4irxl0 


-16 


5x10 


-10 


e.S.U. 


or on in commonly employed units: 


E « 2 volts (2.24) 

P cm (mm Hg) 

This is called low field case. For high field the limit becomes: 


E >5 . 2 volts (2.25) 

p cm (mmHg ) 

The cases in which limit (2.24) holds the mobility of ions are 
independent of electric field so that the drift velocity: 


U, = kE 

a 


To calculate the drift velocity, it is considered that during the 


time between collision, t-^, the ions move a distance d in the field 
experience a force F = eE and an acceleration a = ^ , therefore. 


and 


, 1 2 eE 2 

d - 7 a " 3T c x 


and: 


U 


d_ _ 1 eE 
t x “ 2 m 


EeA 

2mc 


where c is given as: 
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1 

( 8KT \ 2 

c ‘\™) 

thus, the mobility becomes: 

k - (k) | < 2 - 2 » 

Equation (2.26) is a very simple expression for mobility, and 
neglects such effects as role of the mass of the parent gas, and forces 
between ion and gas molecules. There have been numerous theories to 
correct these deficiencies. A thorough discussion can be found in 
Loeb (1955). Langevin derived an expression, based on kinetic theory, 
for mobility as: 


k = 0.815 ~ (2.27) 
Me V m 

where c is the root-mean- square-velocity of the gas molecules . 

A more general treatment which also include high field cases was 
given by Wannier (1953) - He assumes that collision cross section is 
velocity dependence and has the form: 

o(u)u a = r (2.28) 

where a and T are constants* This assumption has two important special 
cases which arise, taking a = 0 and a « 1* The case of a = 0 is that 
of a constant mean free path representing solid elastic impact. The 
case of a - 1 is that of a constant mean free time and is applicable 
to the polarization force. Based on the above assumption, using the 
Boltzmann equation he shows that for low electric field the mean drift 


velocity can be written as: 
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1 

/ \ 2-a 

U d = const j (2.29! 

where N is the number density of molecules . This gives the variation 
of the drift velocity with the field E. 

For constant mean free path, a = 0 (F » a) and equation (2.27) 
reduces to: 




i i 
2 1 2 

= const a A 



(2.30) 


where the drift velocity is proportion to 
to For constant mean free time, a = 

Je 


/E and mobility proportional 

1 (r.= au) : 


/ U \ = const aT = const — T (2.31) 

\ d/ m 


where the drift velocity is proportional to E and mobility is constant. 

Note that in deriving equation (2.30) and (2.31), the quantities 

A = and T = — have been considered. 

No NOU 

For high electric field, the mean drift velocity is given by 
Wannier as : 


a-1 


eE 


( U d> ‘ COnSt miff 
with the special cases: 


(?) 


(2.32) 


1 

2 


< u d) * const a x (ir) 

for constant mean path 0 ~ 0; and: 


„ eE , 
= const — A 
M 


(?) 


1 

2 


(2.33) 
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< U = const a t = const x 

d/ m 

for constant mean free time of a = 1. 

As shown in the above discussion for a = 0 in both high and low 

1 

? 

electric fields the drift velocity is proportional to E and mobility 

1 

-2 

proportional to E . When a = 1 in both cases, the drift velocity is 

proportional to E and mobility is constant, 

Wannier's theory has been successful in interpreting the mobility 

data for many cases. At low applied field the mobility is constant. 

This can be explained as a constant mean free time between ion and 

molecular collisions. Quite frequently, mobility is constant even 

when the limit (2,25) prevails. The region where mobility depends on 
E 

— varies with ion and gases, it probably depends on the strength of 

forces acting between ions and molecules. At still higher field the 

1 

-2 

mobility is proportional to E consistent with a constant mean free 
path between ion and molecular collisions. 

E —1 

For the Gerdien condenser under the investigation — > 2Vcm 

(mm Hg) ^ for an applied potential of 1 volt at altitudes greater than 

52 km. As shown in Figure 2 the electric field for a given potential 

E 

varies within the condenser. It is clear, therefore; that — varies 
from low to high field within the instrument. However, the fact that 
the applied sweep potential changes the electric field from low to high 
field is not important. The important fact is that the mobility (the 
slope of current “voltage characteristic curve) in the condenser remains 


constant. 
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2.5 Fringing Field Effect 

The electric field strength and configuration at the front opening 
of the condenser have a controlling effect on incoming ions or electrons. 
The collecting area of the condenser depends on these stray fields. It 
has been shown (Schmeer et al., 1962) that for large positive bias 
voltages the saturation current, instead of remaining constant, 
actually falls as shown in Figure 3. This happens because lines of 
force originating from the surface of outer cylinder have a repelling 
effect on arriving positive ions. 

It has been suggested (Jespersen et al. , 1968) that for reducing 
this effect both ends of the condenser may be closed by a wire mesh. 

By doing so, the field lines will not extend outside the wire mesh. 

Rose and Widdel (1972), however, have reported this reduces the 
collection efficiency of the instrument considerably (more than 50 
percent) . 

2.6 Photo Emission Effect 

Photo-electrons will be released from metal surface of the probe 
when it is irradiated by solar ultraviolet light and solar x-rays 
(Pederson, 1964) . If sunlight hits the collector of the condenser 
some of the electrons released may be drawn out of the condenser by 
the air flow. This may give rise to a spurious output current, which 
sometimes may be comparable to actual ion current. Shielding can 
prevent this photo emission effect. Photo emission, however, can 
provide the effect of a low impedance connection between an exposed 
electrode and ambient medium. 
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Figure 3. Apparent Decrease of Ion Current as a Result of 
Fringing Field at the Condenser Entrance (after 
Schmeer et al., 1962. 



CHAPTER III 


DESIGN AND INSTRUMENTATION 

3.1 Design of Simplified Version 

A picture of the complete probe is shown in Figure 4. It consists 

basically of two parts: 

(a) An aspiration condenser to measure the ionospheric charge 
particles concentration and their respective mobilities. 

(b) A "blunt probe" assembly which is used to provide the various 
parts of the electronics and telemetry system. 

A schematic cross-section of the whole system is shown in Figure 5. 

The Gerdien condenser is made from two concentric aluminum cylinders. 

The outer cylinder is 10.16 cm long and has a diameter of 7.4 cm and is 
held to the blunt probe by three supporting posts. The inner electrode 
(collector) has a length of half the outer cylinder (5.08 cm) and a 
diameter of 1.4 cm; it is positioned by the center post. This reduction 
in length will prevent the direct sunlight on the collector, therefore, 
it will reduce photo-emission effects. In addition, it reduces the 
fringing effect at the front aperture. Three side posts and the center 
post are electrically isolated with teflon insulation from the outer 
and inner electrodes, (see Figure 5) and they have the same potential 
as the guard ring of the blunt probe. The outer cylinder which is used 
as the driving electrode is electrically connected to the return elec- 
trode of the blunt probe. The collector is connected to the electro- 
meter and has the same potential as the collector disc of the blunt 
probe, therefore, the supporting posts and guard ring are basically of 
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Figure 4. Gerdien Condenser/Blunt Probe Assembly. 

ORIGINAL PAGE IS 
OF POOR QUALITY 






the same potential as the collector in this configuration. 

In order to keep the wiring at a minimum and avoid r-f inter- 
ference, the electronics for the electrometer are seated directly 
behind the guard ring and are shielded in a small aluminum cup. The 
remaining electronic subsystems and telemetry transmitter are contained 
inside the antenna can. Erie RFI filters are used for electrical 
connection through the can to the electrometer. A nylon spacer 
electrically separates the collector from the larger return electrode, 
which is comprised of the aluminum antenna can and power supply wrapped 
with copper tape (Figure 4) . The antenna can is mounted to the power 
supply by cylindrical spacers. The other side is attached to the base 
plate which connects the payload to the parachute and provides the 
mating of the payload with the nose cone and rocket motor. 

This Gerdien condenser /blunt probe assembly, which has a height 
from base plate to the top of outer cylinder (see Figure 4) of 39.37 
cm and weighs 4.5 pounds, is compatible for flying on the Areas rocket, 
a vehicle used by the Meteorological Rocket Network (MRN) . 

3 . 2 Instrumentation 

Figure 6 shows a block diagram of the Gerdien condenser system. 

A brief discussion follows of the sub-system diagrams shown in Figures 
7 to 10. 


3.2.1 Sweep Voltage Generator 

By applying a voltage waveform across the two electrodes of the 
Gerdien condenser and measuring the current collected, it is possible to 
calculate the charge particle concentrations and mobilities in the 
surrounding medium. The voltage waveform used for the Gerdien condenser 
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Figure 8. Logarithmic Electrometer Circuit 
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in this particular experiment and electronics used to generate it is 
shown in Figure 7. 

The sweep voltage generator is basically a free running (astable) 
multivibrator followed by a constant current run-down circuit which 
provides a negative-going ramp function for the output voltage when 
the output of the multivibrator is negative. The output of the multi- 
vibrator has a positive on-time determined by with an amplititude 

of +5 volts followed by a negative pulse of amplitude -5 volts whose 
time duration is controlled by R^. Tlie transistors in the circuit 
are contained on a single integrated circuit chip (RCA3018A) . 

3.2.2 Electrometer 

The electrometer is a highly sensitive ammeter used to measure 
very small currents to the collector of the instrument. A circuit 
diagram and its output response is given in Figure 8. 

The electrometer operational amplifier is an Intersil hybrid 

(ICH8500A), which has extremely low input bias current. The dual 

diodes and D 2 (FJT2000) generate the logarithm response for positive 

10 

and negative output voltages, respectively. The resistor Rp(10 ft) 
provides a linear characteristic for small currents. The 1N459A 
diodes are used to protect the MOSFET input of the electrometer. 

3.2.3 Voltage Controlled Oscillator (VCO) 

The voltage controlled oscillator, shown in Figure 9, transforms 
the output voltage of the electrometer into a pulse train waveform 
suitable for modulating the 1680 MHZ transmitter tube. A simple PPM-AM 
modulating system is used in telemetering the data back to ground. 

The operational amplifier at the input is used as a level shifter. 
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The output of this circuit is pulses with constant amplitude of -5 volts 
whose on-time is determined by R 2 C 2* The de P end s on the 

charging rate of which is linearly proportional to the collector 
current of transistor T^. The use of emitter follower (transistor T 2 ) 
at the output is to reduce the off time and have low output impedance. 

3.2.4 Telemetry System 

The telemetry system makes use of the RCA 1680 MHZ integral tube 
and cavity transmitter which directly feeds a 0.925A. circular slot 
antenna. The voltage controlled oscillator pulses off the transmitter 
at a rate of 0-200PPS. This form of modulation is compatible with the 
GMD/TMQ-5 receiving system, designed for radiosonde balloon data 
acquisition and is commonly used by the Meteorological Rocket Network 
(MRN) . A circuit diagram of the transmitter is given in Figure 10. 

A detailed discussion of the antenna characteristics is given in 
Ionosphere Research Laboratory Scientific Report No. 249 (Cuff in, 

1965). 


3.2.5 Power Supply 

Four 1.25 volts series-connected nickel cadmium batteries provide 
voltage for operating the previously mentioned Gerdien condenser 
electronics- excluding the transmitter tube and cavity. A (DC to DC) 
converter is used to provide +130 volts and -20 volts for the 


electronics . 
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CHAPTER IV 

AERODYNAMICS CONSIDERATIONS 


4.1 Flow Regimes 

The classical theory of the Gerdien condenser discussed in Chapter 
2 is based on the assumption that the instrument operates in a continu- 
ous medium. The continuous regime obtains when the mean free path is 
much smaller than the characteristic length which for the Gerdien 
condenser is the diameter of the outer cylinder. Stated more precisely, 
the collision-controlled regime occurs when the Knudsen number K. <_ .01, 

where K is the ratio of mean free path to the diameter Conley (1974). 
n 

The mean free path calculated (Maris, 1929) is plotted against the 
altitude in Figure 11. Figure 12 shows the flow regime for the Gerdien 
condenser under the investigation. Below 67 km, as is shown in the 
Figure, the mean free path is small compared to the diameter of the 
condenser, and therefore gas dynamics is sufficient for the aerodynamic 
calculation. Since the experimental data obtained in this work lies in 
this height region, the aerodynamic property of the condenser operating 
in this altitude range will be discussed in detail in the next sections. 
However, it is appropriate here to describe briefly the fluid dynamic 
characteristics of the other regimes shown in Figure 12. 

In the region between about 67 km and 102 km the mean free path is 
small but not neglegible compared to the characteristic length of the 
condenser. The gas molecules velocity, in contrast to the continuous 
regime, is no longer zero on the solid surface of the instrument. The 
gas molecules either move on it with a definite velocity or are reflected 
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Figure 11. Variation of Mean Free Path with Altitude 
(after Maris, 1929). 
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Figure 12* The Flow Regimes for the Gerdien Condenser, 



from the surface or both, depending upon the smoothness of the solid 
surface (Kennard, 1938). An ordinary treatment of the problem is no 
longer satisfactory, the problem must be solved based on the kinetic 
theory of the gases . 

In the free molecular regime (above 124 km) the mean free path is 
much greater than the characteristic length of the condenser. An 
analysis based on a Maxwellian velocity distribution for this situation 

seems satisfactory (Tsien, 1946). 

The mean free path in the transition flow is of the same order of 
magnitude as the diameter of the outer cylinder, therefore, the . 
collisions between molecules are of comparable importance to the colli- 
sions of molecules with the surface. The problem is extremely 
complicated and a satisfactory theoretical solution has not yet been 
obtained. 

Therefore, when the probe is operating in a region where the 
continuous flow concept is not applicable, a theoretical calculation 
of velocity distribution in the condenser is very complicated. Although 
by using a condenser with larger diameter, it is possible to extend the 
altitude region where the continuous flow concept is applicable, there 
are practical limitations to this. 

4.2 Compressibility Effect 

As pointed out in the last section the gas dynamic treatment of the 
flow in the collision— controlled regime is sufficient to explain the 
flow properties of the condenser. For air with subsonic velocity and 
small density variation with pressure and temperature, the compressibil- 
ity effect is usually neglected. To a first approximation the 
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compressibility error can be estimated with the relation: 


n 

o po 
n ~ p 

which results from a consideration of the density change in the medium 
caused by adiabatically bringing it to rest (Liepmann and Roshko, 1957). 
In this equation n is the number density; p is the mass density; the 
subscript zero denotes the condition far from the probe; 3 = 1,4 is a 
constant; and M is the Mach number which is the ratio of the velocity 
of the probe-parachute assembly to the velocity of the sound . The 
compressibility error is very small at lower altitudes (about 0.09 
percent at 30KM) and increases with altitude (about 12.5 percent at 
67 km). 

Whether the air flow through the condenser is laminar or turbulent 

depends on the Reynolds number. It can be derived from the general 

definition of Reynolds number, — , where h is hydraulic mean depth 

(defined as the area of the cross-section divided by the perimeter); U 

y 

is the average velocity over a cross-section; and V - — is the kinematic 

K 

viscosity. Thus the Reynolds number for the Gerdien condenser isj 



r e = 


2 (R -R.)U 
Q r 

V 


(4.1) 


Using the probe descent velocity (U) from Figure. 13 and the values of 
coefficient of viscosity (y) and mass density (p) from U, S„ 

Standard Atmosphere Suppliments (1966), the Reynolds number R has been 
plotted against altitude in Figure 14 , If the Reynolds number is less 
than a critical Reynolds number, the flow is laminar. The critical . 
Reynolds number must be determined experimentally; its value for two 






3 
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Figure 14. Variation of the Gerdien Condenser Reynolds 
Number with Altitude. 
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•yoxy long concentrrc cylinders ts 2000 (Bird, Stewart and Lightfoot, 

1960). If this number be adopted as the critical Reynolds number for 
the Gerdien condenser, then the air flow through the condenser is 
laminar (see Figure 14). 

4.3 Velocity Distribution 

The flow condition in a short Gerdien condenser (L = 10.16 cm) 
corresponds roughly to that in the "inlet length" of a long pipe and 
is a highly complicated function of Reynolds number and tube length 
(Paltridge, 1965). Figure 15. a shows schematically the flow condition 
in the inlet of a long pipe. At the entrance a very thin boundary 
layer is formed around the wall of the pipe. The thickness of this 
layer increases downstream in the pipe until it becomes equal to the 
radius of the pipe. Until this happens there is a core of fluid 
unaffected by viscosity. Since the boundary layer thickness increases 
as air flows deeper into the pipe, and because of constancy of flux at 
any cross-section, this core is accelerated. When the viscosity effect 
occupies the whole cross-section of the pipe, a fully established 
laminar flow will be developed which has a parabiolic velocity distribu- 
tion (Shames, 1962). By analogy, the flow condition at the inlet of 
two long concentric cylinders has been shown in Figure 15. b. Note that 
in case of the concentric cylinders two different boundary layers are 
present; however, as long as there exists a core uninfluenced by 
viscosity, the velocity distribution in the core and boundary of outer 
cylinder must be the same as of a pipe with similar size and geometry. 

To determine the inlet length (core length) of flow in the 
condenser the so called "Schiller's Method" will be followed. A detailed 
discussion of this method can be found in Goldstein (1938) . Assuming 




Figure 15. Flow Conditions: (A) in Inlet of a Pipe, (B) in 
Inlet of Concentric Cylinders (after Shames, 
1962). 
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that boundary layer has a thickness 6 at any cross-section the relation: 



can be written. Where is the velocity in the core and y = R q - r . 
This relation along with equation of continuity and momentum equation 
give the following relation between x and U c : 

x 

RR =f (X) 
o E 

and : 

U 

X = - 1 (4 . 2) 

U 

P 

where x is the distance measured from the entrance; is the core 
velocity; and is the average velocity in the fully stabilized 
laminar flow region. Figure 16 shows the variation of f (X) with X 
which has been reproduced from Goldstein (1938). In case of flow in a 
pipe it is well known that the maximum velocity U max is twice the 
average velocity in the parabolic region; then by letting = T ~ T max 

p 

= 2U , y = 1 and the value of f (x) can read from the figure (f (x) * 

P 

.0575), and hence: 

x = .0575 R R 
o E 

In case of concentric cylinders the equality U^ x = 2U^ is no 

longer true, however, a relation between U and U can be found 
07 max p 

assuming a steady incompressible laminar flow through the condenser. 

The Navier-S tokes equation in cylindrical coordinates can be written as 
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iV i??* 

9r 2 r 8r 


1 dP 
U dx 


(4.3) 


which integrates to: 


3r dx 2p r 


(4.4) 


integrating again: 


HP r 

u = - T~ 777 + c i ltl r + S 
x dx 4h 1 2 


(4.5) 


where c and c 2 are constants of integration and can be evaluated by 

applying boundary conditions U = 0, at r = R and at r - R : 

X u i 

R 2 
dP o 

c, InR + = — 7 - 

1 o 2 dx 4p 

R 2 
dP i 

c n In R. + c = r- 77 
1 l 2 dx 4U 


or; 


2 2 
R o- R i dP 
R dx 
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W in 
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(4,6) 
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R 2 In R - R 2 In R. 
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dx 


(4.7) 


The average velocity in the parabolic region is: 

2tf 


it. 


rdrd© 


u = u = 

P X 


2'fT R 

1 1 


rdrdG 


(4.8) 
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which after substituting the value of U from equation (4.5) with c. 

X -L 

and as given in equations (4.6) and (4.7) integrates to: 


U = U = 
P x 


9 2 2 2 2 2 

R + R7 R - R7 / R In R - RT’ In R. . 

o i . o i l oo i l 1 

8(1 4 „ , M ^ 2 - r - " 2 

ln r7 \ 01 


R 2 ln R - R 2 ln R. 

X o o X 


4U ln^ 

X 


The maximum velocity U max can be found by solving 
substituting this value of r in equation (4.5): 


= 0 for r, and 


2 2 

R - R7 
o x 1 

R R 

. o max 

4 U In — 

1 


max dF 
2U dx 


2 2 
R - R7 
o x 

R 

2 ln R^ 

X 


-R 2 ln ~ + ( r 2 - R 2 \ In R + R? ln R - 
max Rj V o x/ max x o 


2 2 
R + R7 
o x 


2 2 
R R - RT 
- o o 1 

ln ir — 2 — 

X 


therefore: 




If rhe value R ^ 3 • 7 cm and R. **• 0*7 cm be inserted, in equation (4.8) , 
o 

then x* 0-54 and f(x) can be read from Figure 16 (f (x) = 0,017). 

With this value of f(x) the equation (4.2) reduces to; 


X = R Q R E f(x) ^0.0629^ (4-10) 

Equation (4.10) approximately determines the distance required for 

fully stabilized laminar flow to be developed. The variation of x 

with altitude has been plotted in Figure 17. Below an altitude of 

about 56 km the inlet length exceeds the length of the condenser, above 

this height range both a core of fluid and a fully stabilized flow 

region exist. Because of nonlinearity in the differential equation, 

only a numerical calculation of velocity in the inlet is possible. 

Goldstein (1938), using the equation of motion and the continuity 

equation has numerically obtained the fluid velocity at every point in 

the inlet of a pipe. Figure 18 has been reproduced from Goldstein; in 

x 

case of the Gerdien condenser the values of velocity for — less than 

o 

about 0.2 must be disregarded. 

It has been suggested (York, 1974) that above 56 km the air 
velocity at the entrance may not be equal to free stream velocity; 
using Bernoulli's equation the velocity at the entrance can be 
approximately determined: 
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Figure 17. Variation of Flow Core-Length with Altitude. 
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+ P 
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‘ atm 
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(4.11) 


where U. and P. are the air velocity and air pressure at the front 
mm 

end of the condenser, respectively; I*p ro ^ e *-" s * : t he instrument velocity 
(free stream velocity) . The pressure at the entrance may be obtained 
from: 


P. 

xn 


= p 

atm 


+ ap vl 


(4.12) 


where Ap. is the pressure loss in the condenser due to the air 
VL 

viscosity and is given by Shames (1962): 


Ap 

VL 2 D P 


(4.13) 


where f is the friction coefficient and for laminar flow is given by: 


f 



(4.14) 


substituting equation (4.12) for P. in equation (4.11) with AP and 

in VL 

f as are given in equations (4,13) and (4.14), respectively, the 
equation (4.12) can be solved for U. : 


U t 

in 
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+ U 
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probe 



(4.15) 


If it can be assumed that the average velocity is the same as the free 

stream velocity, U, can be calculated using an atmospheric model for 

in 

pressure and density (CIRA) . It turns out, however, there is no 
significant difference between the values of IL calculated from equation 
(4.15) and free stream velocity. From the above discussion, it can be 
concluded that in this case (above 56 km) the air enters the condenser 
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with a velocity of about the same as free stream velocity. Within the 
condenser two distinct flow regions exist: 

(a) An entry region with a relatively thick boundary layer in 
which the air flow velocity is a complex function of x and 
r; and a core of fluid with air velocity as a function of x 
only. 

(b) A fully stabilized flow region with a parabolic velocity 
distribution which its average is equal to the free stream 
velocity. 

Since an overall average velocity calculation is not possible, it is 
arbitrarily assumed that average velocity at any cross-section is the 
same as the free stream velocity. 

When the boundary layer thickness is thin (below 56 km) , the so 
called "displacement thickness" concept may be applied. For the outer 
cylinder the boundary layer thickness is much smaller than the radius, 
and a flat plate approximation may be used: 



(4.16) 


For the inner cylinder the displacement thickness is given by Kelly 


(1954) : 



Since 6 and <5. are a function of x, their average over the entire 
o x 

length of the condenser may be obtained as: 
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(4.19) 


Using the constancy of flux concept at any cross-section the core 
velocity can be calculated: 


tt[(R - 5 ) 2 - (R. + 6 ) 2 ]U = TT (R 2 - R 2 )U (4.20) 

1 v >o o x i c o X o 

A simple iterative computer program given in Appendix A has been used 

to determine the values of U as a function of altitude; the results 

c 

of this computer calculation is plotted in Figure 19 (dotted line) * It 
must be noted, although; the average velocity in the condenser is 
increased considerably with comparison to the descent velocity, it will 
not alter the density and mobility of charge particles using the free 
stream velocity (since the product of velocity and area remains 
constant in this method) . 



- 50 - 



Figure 19 . Descent Velocity of the Gerdien Condenser (full line) 
and Calculated Velocity of Air Flov (clotted line) . 



CHAPTER V 


FLIGHTS AND DATA 

5.1 The Gerdien Condenser Launch 

The Gerdien condenser and associated instrumentations, which were 
described in chapter three, is compatible for flying on the Areas 
rocket, a vehicle routinely used for wind and temperature measurement. 
This rocket is designed to carry a parachute-borne payload for ejection 
near Apogee. Three parachute-borne Gerdien condenser experiments 
employing the Areas rocket as the launch vehicle were carried out. 

Most of the useful data obtained from these experiments occurs after 
expulsion from the rocket when the Gerdien condenser is descending 
subsonically to the earth suspended under a parachute. 

The first successful Gerdien condenser flight occurred at White 
Sand Missile Range (WSMR) , New Mexico on 11 January 1974 at 1845 UT. 

The flight apogee for this experiment was about 70 km. An analysis of 
data indicates that the angle of attack problems have been largely 
eliminated by the use of a more aerodynamic ally stable decelerating 
device, a starute manufactured by Space Data Corporation. 

The other two flights took place at Wallops Island (WI) , Virginia; 
one on 29 June 1974 at 2112 UT and the other on 30 June 1974 at 0315 
UT. The data from these last two experiments have not yet been reduced 
and will not be presented in this work. 

The probe is activated and calibrated on the ground just prior to 
launch; the preflight calibration procedure involves inserting a 
precision resistor between two electrodes of the probe and recording 
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the telemetered probe response on a magnetic tape* The calibration 
response of the Gerdien condenser to a sweep voltage waveform (Figure 
21. a) is shown in Figure 20. While the probe is in— flight, the impedance 
of the surrounding medium (or conductivity) is the unknown quantity 
which is measured; then knowing saturation current and volume rate of 
air flow through the aspiration condenser both density and mobility can 
be determined (see Chapter II). The in-flight mode of operation (a 
typical strip chart record) at an altitude of about 43 km is shown in 
Figure 21* bo 


5*2 Reduction of Data 

The first step for reducing data is to draw the current-voltage 
characteristics from the strip chart record which represents the 
variation of frequency with time* This involves obtaining current 
from calibration response (Figure 20) and voltage from sweep voltage 
waveform (Figure 21. a) for different points on the strip chart* Once 
the current-voltage characteristic is obtained, a radar plot showing 
the position of the payload as a function of time (Figure 22) provide 
the information necessary for determining the altitude at the time of 
the measurement* 

When sweeping the voltage to outer cylinder in the Gerdien 
condenser, a displacement current will go to the electrometer. This 
current is given by the equation 


i 


d 



(5.1) 


where c is capacity between central electrode and electrometer and 
is constant for a ramp voltage sweep. The displacement current was 
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Figure 20. Calibration Current as a Function of Frequency. 
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Figure 21. The Gerdien Condenser Sweep Waveforms: 
(a) Applied Sweep Voltage, (b) Strip 
Chart Record . 
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calibrated before flight and has a value of 1.2 PA. 

Figure 23 shows two different current-voltage characteristics of 
the Gerdien condenser obtained from the strip chart record for altitudes 
of about 62 km (Figure 23. a) and 39 km (Figure 23. b). The effect of 
displacement current has been taken into account in drawing the current- 
voltage characteristic; it has been added to a positive collector 
current and has been subtracted from a negative one. Two distinct 
slopes of the current-voltage characteristic in Figure 23. a indicate 
the presence of two groups of ions with different mobilities; the 
current-voltage characteristic of the Gerdien condenser above about 47 
km is of this type. The Gerdien condenser characteristic in altitude 
range between 38 km and 47 km is similar to the one shown in Figure 
23. b; hence, a spectrum of positive ions with different mobilities 
must be present in this altitude range. 

For calculating the positive ion densities and their relevant 
mobilities, the equations: 


and : 


N = 


.7 2 - 

elKR" - R,)U 
o i 


for V > V 


di. 


dv: 


k = ~r — “ for V < V 

N e — 
e 


(5.2) 


(5.3) 


haye been used. The equation (5.2) is obtained from equation (2.-11) 
with V q as is given by equation (2.12), and equation (5.3) is derived 
from equation (2.10) and (5.2). 

The results of positive ion densities and positive ion mobilities 
as a function of altitude for flight, made on January 11, 1974 at WSMR 







- 58 - 


are given in Figures 24 and 25, respectively. Above 47 km two groups 
of positive ions are present; the heavier group is about an order of 
magnitude less mobile than the lighter group. Between about 38 and 47 
km the positive ions with different densities and their mobility 
spectrum are represented by straight lines. Below about 38 km the 
applied sweep voltage is not high enough to saturate the Gerdien 
condenser, consequently; no data is given. A preliminary analysis of 
the strip chart data indicates that electrons have not been collected 
for the most of the flight, therefore, the negative ion densities and 
mobilities have not been calculated. 
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Figure 25. Positive Ion Mobility as a Function of Altitude. 


CHAPTER VI 


CONCLUSIONS 

The absence of reliable D-region particle mobility data has 
resulted in the design of a simple Gerdien condenser (aspiration 
probe) system believed capable of accurately measuring ion density and 

mobility. The angle of attack problem associated with earlier instru- 

* 

ment has been largely reduced by the use of an aerodynamically more 
stable device, a starute, manufactured by Space Data Corporation. 

Because of nonlinearity in the differential equations involved, it 
has not been possible to determine the average velocity of air flow 
through the condenser theoretically, even in case of relatively simple 
incompressible laminar flow (Chapter IV). Therefore, the average veloci- 
ty or the volume rate of air flow through the Gerdien condenser must be 
evaluated experimentally; it can be measured using a flow-meter (such 
as the one used by Rose and Widdel, 1972). 

Experimental results of 11 January 1974 at WSMR (Figures 24 and 
25) show the positive ion densities and mobilities between an altitude 
range of about 38 and 65 km. At lower altitude the applied sweep 
voltage was not high enough to saturate the condenser; it seems 
increase of positive peak of voltage waveform from +5 to +10 volts is 
sufficient to saturate the instrument having the descent velocity of 
this order (Figure 13). Electrons were not collected effectively in 
this experiment. A possible reason may be electrons having relatively 
high mobility were reflected before reaching the device. 
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APPENDIX A 
COMPUTER PROGRAM 


A SIMPLE PROGRAM TO CALCULATE DISPLACEMENT THICKNESS 
IN 'GERDIEN CONDENSER’. 

A=3 « 7E-2 
B=7E-3 
ABS-1E-3 
WRITE (6, 16) 

GO TO 20 

20 READ(5,10,END=100)M,V,U0,D 
10 F0RMAT(12,E10,3, 2X,F10. 1, 5X,E8. 2) 

UC-UO 
GO TO 14 

14 P=SWRT(V*D/UC) 

DEL0=3, 442/3. *P 

'DELO' IS THE AVERAGE DIS. THICKNESS IN OUTER CYLINDER, AFLAT PLATE 
APPROX, HAS BEEN ASSUMED. 

R=4./B*P 

DELI-1. 721*P* (2. /3 . 102*R+. 292/5 .*R**2) 

’DELI’ IS THE AVERAGE DIS. THICKNESS ON INNER CYLINDER, BASED 
ON KELLY'S RESULT, 

F=3 , 14* (A**2-B**2) *UO 

FDEL=3 o 14* ( (A-DELO) ** 2- (B+DELI) **2) *UC 

F AND FDEL ARE FLUXES AT ENTRANCE AND ACROSS SECTION, WHERE 
DIS. THICKNESS HAS BEEN TAKEN INTO ACCOUNT, RESPECTIVELY. 
IF(FDEL.GT.F) GO TO 12 
UC=UC+ABS 

A ITERATIVE METHOD IS USED, UC IS INCREASED BY AN INCREMENT 
OF .001 UNTIL FDEL-F (CONSTANCY OF FLUXES). 

GO TO 14 

:.12. WRITE (6, 18)M, DELI, DELO, UC,F, FDEL 
18 FORMAT (15X,12,2(8X,E10. 4) ,10X,F6. 2, 2 (9X,E10 . 4) ) 

GO TO 20 

16 FORMAT (13X, 'ALT, : ’ ,13X, 'DELI: ’ ,13X, 'DELO: ' ,6X, 'CORE VEL.:’,5X, 
l’FLUX AT ENT.: ',8X,' CORE FLUX:’//) 

100 STOP 
END 
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